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Abstract 

Characteristic examples are presented of scenarios of particle production and decay 
in supersymmetry models in which the supersymmetry breaking is transmitted to the 
observable world via gauge interactions. The cases are chosen to illustrate the main 
classes of LHC phenomenology that can arise in these models. A new technique is 
illustrated that allows the full reconstruction of supersymmetry events despite the 
presence of two unobserved particles. This technique enables superparticle masses to 
be measured directly rather than being inferred from kinematic distributions. It is 
demonstrated that the LHC is capable of making sufficient measurements so as to 
severely over-constrain the model and determine the parameters with great precision. 



*This work was supported in part by the Director, Office of Energy Research, Office of High Energy 
and Nuclear Physics, Division of High Energy Physics of the U.S. Department of Energy under Contracts 
DE-AC03-76SF00098 and DE-AC02-98CH10886. 



Contents 

1 Introduction 



2 Point Gla 

2.1 Lepton/Photon Distributions 

2.2 Reconstruction of G Momenta 

2.3 Reconstruction of Gluinos and Squarks 

3 Point GT5 

3.1 Effective Mass Analysis at Point GTE . 



3.2 Selection of Dilepton events" 



3.3 Extraction of fluino 


3.4 Extraction of squarks 




4 Point (42a 




4.1 Dilepton Distributions 


4.2 Detection of on 


4.3 Detection of xt decay 


3 



5 Point G2b 

|5.1 Effective mass analysis at Point G2b . . . 

|5.2 Slepton mass determination 

|5.3 Reconstruction of y 1 /, y!j and~\^ 

[5.4 Extraction of Tj, 

5.5 Reconstruction of Squarks 

5.6 Reconstruction of tt Decays 

5.7 Direct production of electro- weak sparticles 



6 Determining SUSY parameters 



7 Discussion and Conclusions 



1 Introduction 



If supersymmetry (SUSY) exists at the electroweak scale, then gluinos and squarks will be 
copiously produced in pairs at the LHC and will decay via cascades involving other SUSY 
particles. The characteristics of these decays and hence of the signals that will be observed 
and the measurements that will be made depend upon the actual SUSY model and in 
particular on the pattern of supersymmetry breaking. Previous, detailed studies of signals for 
SUSY at the LHC |], §, || have used the SUGRA model [§, |5| , in which the supersymmetry 
breaking is transmitted to the sector of the theory containing the Standard Model particles 
and their superpartners via gravitational interactions. The minimal version of this model 
has just four parameters plus a sign. The lightest supersymmetric particle (x?) has a mass 
of order 100 GeV, is stable, is produced in the decay of every other supersymmetric particle 
and is neutral and therefore escapes the detector. The strong production cross sections and 
the characteristic signals of events with multiple jets plus missing energy fa or with like- 
sign dileptons £ ± £ ± plus fa @ enable SUSY to be extracted trivially from Standard Model 
backgrounds. Characteristic signals were identified that can be exploited to determine, with 
great precision, the fundamental parameters of these model over the whole of its parameter 
space. Variants of this model where R-Parity is broken |7| and the Xi is unstable have also 
been discussed |§. 

There also exists a class of gauge mediated SUSY breaking (GMSB) models P, |TU 
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which the supersymmetry breaking is mediated by gauge interactions. The model assumes 
that supersymmetry is broken with a scale yF in a sector of the theory which contains heavy 
non-Standard-Model particles. This sector then couples to a set of particles with Standard 
Model interactions, called messengers, which have a mass of order M. These messengers 
are taken to be complete representations of SU(5) so as to preserve the coupling constant 
unification of the Minimal Supersymmetric Standard Model (MSSM). The mass splitting 
between the superpartners in the messenger multiplets is controlled by \f~F. One (two) loop 
graphs involving these messenger fields, then give mass to superpartners of the gauge bosons 
(quarks and leptons) of the Standard Model. This model is preferred by some because the 
superpartners of the Standard Model particles get their masses via gauge interactions, so 
there are no flavor changing neutral currents, which can be problematic in the SUGRA 
models. 

The characteristic spectra of superparticles are different from those in the SUGRA model; 
in particular, the lightest supersymmetric particle is now the gravitino (G). This particle has 
feeble interactions and can be produced with significant rate only in the decays of particles 
which have no other decay channels. In the SUGRA models, G has a mass of order 1 TeV 
and is phenomenologically irrelevant (except possibly for cosmology). In GMSB models, the 
next lightest supersymmetric particle, which we will refer to as the NLSP, decays into G. 
The lifetime of the NLSP is very model dependent: 0(l/xm) < cr < 0(many km). As it is 
not stable, it can either be charged or neutral. 

If the NLSP is neutral, it is the lightest combination (xi) °f gauginos and Higgsinos, and 
it behaves, except for its decay, in the same manner as the x\ in the SUGRA model. (The 
rather unlikely possibility that it is a gluino JiTj] is not explored in this paper.) If its lifetime 
is very long so that none of the produced Xi decay within the detector, the phenomenology 
is very similar to that in the SUGRA models. 
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If the NLSP is charged, it is most likely to be the partner of a right handed lepton. Two 
cases are distinguished. If tan (3, the ratio of the vacuum expectation values of the two Higgs 
fields, is small, then e R , fin and tr are degenerate. If tan (3 is large, f R is the lightest slepton 
and the others can decay into it with short lifetimes. If the lifetime of the NLSP is very 
long, each SUSY event contains two apparently stable charged particles [Q. If it is short, 
each event contains two charged leptons from each of the decays. 

The simulation in this paper is based on the implementation of the minimal GMSB model 



in ISAJET fl3| . The model is characterized by A = F/M m , the SUSY breaking mass scale; 
M m , the messenger mass; N 5 , the number of equivalent 5 + 5 messenger fields; tan/5, the 
usual ratio of vacuum expectation values of the Higgs fields that couple to the charge 2/3 
and 1/3 quarks; sgn/i = ±1, the sign of the \i term, the value being determined by the 
Z mass from usual radiative electroweak symmetry breaking; and C grav > 1, the scale factor 
for the gravitino mass which determines the NLSP lifetime (r NLSP ~ C 2 rav ). At the scale 
M m , for example, the masses of the gluino, squark and slepton are given by 

m~ g = ^AA 5 

ml = ^|A 2 iV 5 + A 2 iV 5 



<-L 
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m 2 = j^l A 2 N 

eR 40tt 2 5 

ml = 4,A 2 iV 5 + Ma 2 A 5 + A 2 N 5 
6tt 2 32tt 2 480tt 2 

ml = ^A 2 A 5 + ^L\ 2 A 5 
UR 6tt 2 5 30tt 2 5 

Here a s , a.i and a.\ are the coupling constants of SU(3), SU(2), and U{\) respectively. 
These masses are then evolved from the scale M m to the weak scale, inducing a logarithmic 
dependence on M m . As in the case of the SUGRA models, this evolution leads to the 
spontaneous breaking of electro-weak symmetry as the large top quark Yukawa coupling of 
the top (and possibly bottom) quark induces a negative mass-squared of the Higgs field(s). 
From these equations it can be seen that as N$ is increased the slepton masses increase more 
slowly than the gaugino masses. Hence only for small values of N$ will the NLSP be a x®, 
at larger values it is a (right) slepton. Messenger fields in other representations of SU(5) 
can be included by modifying the value of N$. If there are several messengers with different 
masses their effects can be approximated by changing N 5 . Therefore, we will consider N 5 to 
be a continuous variable when we estimate how well it can be measured at LHC. 

This paper presents a series of case studies for this model which illustrate its characteristic 
features. We use the four sets of parameters shown in Table 1. These cases are paired and 
differ only in the value of C grav and hence in the lifetime of the NLSP. The masses of the 
superpartners of the Standard Model fields are given in Table 2. In cases Gla and Gib, the 
NLSP is X\] its lifetime is quite short, cr = 1.2 mm, in the former case and long, cr = 1 km, 
in the latter. In cases G2a and G2b the NLSP is a stau. In the latter case it has a very 
long lifetime, cr w 1.25 km, and exits a detector without decaying. The decay £ — > t\t£ 
is not kinematically allowed, so the cr and jln are also long lived. In the case of point 
G2a, the fx, en and jiR are short lived with cr = 52 fim. The production cross section 
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Table 1: Parameters for the four case studies in this paper. 



Point A M m N 5 tan/5 sgn/i C grav > 1 
(TeV) (TeV) 



Gla 


90 


500 


1 


5.0 


+ 


1.0 


Gib 


90 


500 


1 


5.0 


+ 


10 3 


G2a 


30 


250 


3 


5.0 


+ 


1.0 


G2b 


30 


250 


3 


5.0 


+ 


5 x 10 3 



Table 2: Masses of the superpartners, in GeV, for the cases to be studied. Note that the first 
and second generation squarks and sleptons are degenerate and so are not listed separately. 



Sparticle Gl G2 Sparticle Gl G2 



9 


747 


713 








xt 


223 


201 


xt 


469 


346 


X? 


119 


116 


xl 


224 


204 


xl 


451 


305 


xl 


470 


348 




986 


672 


UR 


942 


649 


d L 


989 


676 


d R 


939 


648 


h 


846 


584 




962 


684 




935 


643 


h 


945 


652 


e L 


326 


204 




164 


103 




317 


189 


T 2 


326 


204 


T\ 


163 


102 


V T 


316 


189 


h° 


110 


107 


H° 


557 


360 


A 


555 


358 


H± 


562 


367 



for supersymmetric particles is quite large; 7.6 pb and 22 pb at lowest order in QCD for 
cases Gl and G2 respectively. Note that the larger value of N$ in the case of G2 results in 
considerably smaller squark masses and hence the larger cross section; the gluino masses are 
very similar. 

All the analyses presented here are based on ISAJET 7.37 [13] and a simple detector 
simulation. At least 50K events were generated for each signal point. The Standard Model 
background samples contained 250K events for each of ti, WZ with W —>■ eu, fiu, ru, and Zj 
with Z — > uv, tt, and 5000K QCD jets (including g, u, d, s, c, and b) divided among five 
bins covering 50 < px < 2400 GeV. Fluctuations on the histograms reflect the generated 
statistics. On many of the plots that we show, very few Standard Model background events 
survive the cuts and the corresponding fluctuations are large, but in all cases we can be 
confident that the signal is much larger than the residual background. The cuts that we 
choose have not been optimized, but rather have been chosen to get background free samples. 

The detector response is parameterized by Gaussian resolutions characteristic of the 
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ATLAS [[14]] detector without any tails. All energy and momenta are measured in GeV. In 
the central region of rapidity we take separate resolutions for the electromagnetic (EMCAL) 
and hadronic (HCAL) calorimeters, while the forward region uses a common resolution: 



EMCAL 10%/V£ © 1%, \rj\ < 3 
HCAL 50%/v^ © 3%, < 3 
FCAL 100%/v 7 ^ © 7%, >3 

A uniform segmentation At] = A0 = 0.1 is used with no transverse shower spreading; this 
is particularly unrealistic for the forward calorimeter. Both ATLAS 0] and CMS have 



finer segmentation over most of the rapidity range. An oversimplified parameterization of 
the muon momentum resolution of the ATLAS detector including a both the inner tracker 
and the muon system measurements is used, viz 



5p T /p T = ^0.016 2 + (0.0011p r ) 2 

In the case of electrons we take a momentum resolution obtained by combining the electro- 
magnetic calorimeter resolution above with a tracking resolution of the form 

0.4 



Spt/pt = I i + (3 q^Tji ) v (°- 0004 ^) 2 + °- 0001 

This provides a slight improvement over the calorimeter alone. Missing transverse energy is 
calculated by taking the magnitude of the vector sum of the transverse energy deposited in 
in the calorimeter cells. 



Jets are found using GETJET |T3[ with a simple fixed-cone algorithm. The jet multi- 



plicity in SUSY events is rather large, so we will use a cone size of 



i? = V /(A^) 2 + (A0) 2 = O.4 

unless otherwise stated. Jets are required to have at least > 25 GeV; more stringent 
cuts are often used. All leptons are required to be isolated and have some minimum p? and 
\t)\ < 2.5. The isolation requirement is that no more than 10 GeV of additional be present 
in a cone of radius R = 0.2 around the lepton to reject leptons from 6-jets and c-jets. In 
addition to these kinematic cuts a lepton (e or //) efficiency of 90% and a 6-tagging efficiency 



of 60% is assumed fi~4]| . Where relevant, we include the possibility that jets could appear as 
photons in the detector due to fragmentation where most of the jet energy is taken up by 
7r°'s. Jets are picked randomly with a probability of 10~ 3 . They are then called photons and 
removed from the list of jets. This is a conservative assumption, ATLAS is expected to have 
a better rejection 

Results are presented for an integrated luminosity of 10fb _1 , corresponding to one year 
of running at 10 33 cm _2 s _1 so pile up has not been included. We will occasionally comment 
on the cases where the full design luminosity of the LHC, i.e. 10 34 cm _2 s _1 , will be needed to 
complete the studies. For many of the histograms shown, a single event can give rise to more 
than one entry due to different possible combinations. When this occurs, all combinations 
are included. 
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The next sections of this paper contain detailed examples of analyses that could be carried 
out in the selected cases. In particular, we illustrate a technique where the momenta of each 
of the G"s can be reconstructed even though only the sum of their transverse momenta is 
measured directly. Using these momenta we are then able to reconstruct the squark and 
gluino decays. We devote considerable space to this technique as it is new and enables 
the masses of superparticles to be directly measured rather than inferred from kinematic 
distributions. This technique should be applicable to other cases. We then estimate how well 
the LHC could determine the parameters of the model and comment on possible ambiguites 
in interpreting the signal. In those cases where the signal is characteristally differant from 
the SUGRA cases, the expected precision on the model parameters is even better than that 
in the SUGRA cases. Finally we comment upon how our results can be used to estimate 
how well the LHC would be able to study other parameter sets of the GMSB model. 

2 Point Gla 

In case Gla, the NLSP is Xi an d decays to G7 with cr = 1.2 mm. The total SUSY cross 
section is 7.6 pb. SUSY events are characterized by two hard isolated photons plus the usual 
jets, leptons, and missing transverse energy fx- The presence of two photons in almost 
every event renders the Standard Model backgrounds negligible. The first evidence for new 
physics in this case will arise from a huge excess of events with two photons and missing Et 
over that expected from the Standard Model. 

In a small fraction (2.0%) of the events, the NLSP will undergo a Dalitz decay to e + e~G. 
The electron and positron can be used to accurately determine the decay vertex and a precise 
measurement of the lifetime made. The systematic limit on the precision from the resolution 
of the vertex system of ATLAS is at the 10/xm level. The precision will therefore be limited 
by statistics of the thousand or so observed Dalitz decays. This measurement is important 
as it provides the only constraint on C grav and hence on the true scale of SUSY breaking. 

The effective mass is defined to be the scalar sum of the missing transverse energy Ifi? 
and the pt's of the four hardest jets, which are required to contain more than one charged 
particle with p T > 1 GeV: 

M cS = $ T + p T ,i+ PT,2 + PT,3 + PT,4 ■ 
This is a good measure of the hardness of the event. Events are selected to have 

• M cS > 400 GeV; 

• E T > O.lMeff. 

Electrons and photons are required to have pt > 20 GeV; muons are required to have 
p T > 5 GeV. At least two photons and two leptons are also required for all the analyses in 
this section. 

2.1 Lepton/Photon Distributions 

We begin a detailed study by first selecting events that have at least two leptons (either 
electrons or muons) and two photons satisfying the cuts described above. In order to cleanly 
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Figure 1: M« distribution for e + e + fi + fi — e ± fi T events with two photons and two isolated 
leptons arising at Point Gla. The Standard Model background is negligible. 

select events arising from the decay chain 

we observe that these leptons are correlated in flavor and hence we take events where a 
pair of the leptons have opposite charge and form the flavor subtracted combination e + e~ + 
yU + /i~ — e ± /i T . The distribution in the mass distribution is shown in Figure [l]. This 
distribution has a sharp edge at 



' M In 



105.1 GeV. 



which arises from the decay chain above. If we do not form the flavor combination, the edge 
is still clearly visible but there is much more combinatorial background. The position of this 
edge can be used to determine this combination of masses with great precision given the 
huge statistical sample. The ultimate limit will arise from systematics in the measurement 
of dilepton mass distribution, which given the large sample of Z — > £ + £~ decays that will be 
available for calibration, we expect to have an uncertainty of order 0.1% or 100 MeV. 

From the same sample of events with two leptons and two photons, the smaller of the 
two £ + £~7 masses is formed and the resulting mass distribution is shown in Figure 0. The 
endpoint is at 

'M 2 -M 2 = 189.7 GeV. 

X2 X'j 

the kinematic limit for x° ~~ > G£ + £~^. Instead of having a sharp edge like Figure [l], this 
distribution vanishes linearly because of four-body phase space. The figure shows as a 
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Figure 2: M« 7 distribution for e + e~ + n + jjT — e ± fi T events with two photons at Point Gla. 
A linear fit from 160 to 190 GeV used to determine the position of the 4-body endpoint is 
also shown as a dashed line. 



dotted line a linear fit to the end of the spectrum. From this fit the precise end point can 
be determined. The smaller statistical sample and the need for a fit imply that the resulting 
uncertainty on the value of the end point is larger. We estimate the precision from Figure |^ 
to be ±500 MeV; the full LHC luminosity of 100 fb _1 should enable this uncertainty to be 
reduced to ±200 MeV. 

The subset of events where one mass is larger and the other smaller than this 

endpoint was then selected. Only the combination with the lower iij mass can come from 
X® decay. The M^± 7 distribution for this combination is shown in Figure ^. Two structures 
are present in this plot. There is a sharp edge at 

y M l.- M x? = 112 - 7GeV 

from the photon and the second ("right") lepton in the decay chain (i —>■ £7), and there is 
a distribution that vanishes linearly at 



< / M 2 — M~ = 152.6 GeV 

from the photon and the first ( "wrong" ) lepton. The former of these has a background from 
the latter. The plot shows a linear fit between 115 and 150 GeV, which can be used to 
determine the position of the second endpoint. We estimate the errors on these quantities to 
be ±200 and ±500 MeV for luminosity of 10 fb _1 where they are limited by statistics. We 
expect that that they will become limited by systematics at ±100 and ±200 MeV at LHC 
design luminosity. 
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Figure 3: M^± 7 distribution for e + e~ + fi + \i~ — e ± fi T events sample with two photons; a 
linear fit from 115 to 150 GeV is used to determine the endpoint for the photon and the 
"wrong" lepton and is shown dashed line. 

These four measurements are sufficient to determine the masses of the particles (%2, ^Ri 
and Xi) i n this decay chain without assuming any model of SUSY breaking. Of course the 
existence of, and rate for, this decay chain are model dependent. So is the interpretation of 
the slepton mass as the mass of the £r. A similar combination of three-body and four-body 
distributions will be useful in other cases for which a decay chain involving three two-body 
steps can be identified|] 

2.2 Reconstruction of G Momenta 

The supersymmetry events each have two unobserved G"s. The sum of their transverse 
momenta is, up to resolution effects and possible missing neutrinos, measured as the two 
components of $t- There appears to be insufficient information to reconstruct their mo- 
menta. However the decay chain 



provides three precisely measured final state particles and three mass constraints using the 
masses determined from the measurements in the previous subsection. A 0C fit for the G 
momentum is then possible assuming that Mq = 0. The solution has a 4-fold ambiguity 
since the leptons cannot be uniquely assigned and there is quadratic ambiguity from the 
solution to the constraints. The details are presented in the appendix. 

^For example, at SUGRA Point 5 considered previously Jl], |[, the decay q~L — > X°9 ~ * ^r^ T 1 ~^ Xi^ + ^~<7 
can be used to determine an £ + £~ edge, an £q edge, and an £ + £~q endpoint. 
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Figure 4: The distribution in \Ap\ / \p\ for reconstructed gravitino momenta. The gravitino 
momenta are reconstructed using the method described in the text. The difference (Ap) 
between the reconstructed (p) and generated momentum is then formed for the combination 
with the lowest x 2 ■ Events are included where the lowest x 2 is l ess than 10 or 1 (dashed 
line) . 



We want to have a x® ~ * £ + £~jG decay on both sides of the event. We therefore select 
events with four leptons and two photons. We require two opposite-sign, same-flavor lepton 
pairs and a unique way of combining these with the photons consistent with x\ decay. Then 
there are a total of 16 solutions. fa can be used to resolve these ambiguities as follows. The 
best solution is selected by summing the momenta of the two G"s and calculating the x 2 f° r 
matching this to the measured value of fa, 



X 




fa ~ Ply - Ply 



We then select the solution that has the lowest x 2 (Xmin)- ^ i s assumed that the resolution 
in fa is determined by the total transverse energy E r p, 



A$ x = A fa = 0.6^ E T + 0.03E T 

We can evaluate the effectiveness of this method by comparing the reconstructed values 
of the G momentum to the best match with the generated values. Figure |] shows the distri- 
bution of the fractional difference between the generated and reconstructed G momentum 



Ap 
p 



-Pi 



\P 



gen | 



for solutions where Xmm < 10; the distribution for Xmin < 1 is very similar. As can be seen 
from the figure which indicates a resolution of order 10% with a long tail, showing that the 



10 



150 



o 100 
> 

CD 

u 

o 

<p 

> 50 



" 500 1000 1500 2000 

M zjj (GeV) 

Figure 5: Invariant mass of x° an d two jets, for events at Point Gla. The peak is due to the 
decay g — > qq\\- Events are required to have two x° momenta reconstructed and at least 
four jets with p T > 75 GeV. 

method works quite well. In approximately 40% of the events that enter this analysis, the 
G momenta reconstructed within 10% of their nominal value. The shape of this distribution 
is dominated by detector resolution on the leptons and photons. It is much narrower if the 
generated momenta are used, and it is significantly wider if the resolution for the electron 
energy is taken from the calorimeter alone; the central tracker is important for soft electrons.^ 

2.3 Reconstruction of Gluinos and Squarks 

Squark production is significant at Point Gla even though squarks are considerably heavier 
than gluinos. The reconstruction technique of the previous section can be extended to allow 
the chain q — > gq — > X2QII to be reconstructed and the gluino and squark masses measured. 
We begin with the events selected in the previous section where two X2 momenta each 
arising from X2 ~^ £ + £~jG have been reconstructed. Jets are are then searched for that have 
Pt > 75 GeV in a cone R = 0.4. We require that at least four such jets be present in the 
event. Each x® * s then combined with two and with three jets; the resulting invariant mass 
distributions are shown in Figures [5] and ^ respectively. The jjxt mass of Figure has a 
peak close to the gluino mass, of 744 GeV, while Figure [| has a broader peak near the squark 
masses, 935-985 GeV. The peaks occur in essentially the same place if the jet cut is raised 
to Pt > 100 GeV and so are not simply reflections of the kinematic cuts. It is important to 

Ht is possible in principle to apply this method to other cases involving three identifiable decays, e.g., the 
qL decay chain for SUGRA Point 5 0, || . Unfortunately, the combination of more combinatorial background 
and relatively poor resolution for jets seems to make it not possible to reconstruct the \i momenta in this 
case. 
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Figure 6: Invariant mass of xi an d three jets, for events at Point Gla. The peak is due to 
the decay q — > qg — > qqqx\- 
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Figure 7: Scatter plot of M(x%jjj) vs - M(x23j) f° r events at Point Gla where two X2 
momenta are reconstructed and there are at least four jets with > 75 GeV. 

emphasize that this technique enables the masses of the squarks and leptons to be measured 
directly rather than being inferred from features in kinematic distributions. 

The peaks can be sharpened up considerably by searching for correlations since, for each 
X2 momentum and set of three jets, there should be one X2jjj peak at the squark mass mass 
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Figure 8: Projection of Figure [7| on the M{j$jj) axis for 800 GeV < M(x° 2 jjj) < 1000 GeV. 
The peak corresponds to gluino decay. The fit shown is a Gaussian over the range 500 to 
900 GeV. 
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Figure 9: Projection of Figure ^ on the M^jjj) ax i s f° r 600 GeV < M(y^33) < 800 GeV. 
The peak corresponds to squark decay. The fit shown is a Gaussian over the range 760 to 
1020 GeV. 



and three X2JJ m ass combinations, one of which is at the gluino mass. The scatter plot of 
all combinations is shown in Figure |7|. While the points show significant scatter, there is, 
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nevertheless, a clear peak in the M(x\jj) projection, made by selecting events in the range 
800 GeV < M(x2jjj) < 1200 GeV, shown in Figure ||. The smooth curve shown in this plot 
is a Gaussian fit over the range 500 to 900 GeV; shows a clear maximum at 699 GeV. If a cut 
around this peak, 600 GeV < MijQjj) < 800 GeV, is made a projection of the scatter plot 
onto the jjjx® ax is made as shown in Figure ||, a somewhat narrower peak at the squark 
mass than that of Figure |] is obtained. The smooth curve is again a Gaussian fit over the 
range 760 to 1020 GeV and has a maximum at about 909 GeV. The statistical errors on 
the locations of the peaks are quite small and the precision of the measurements is likely 
dominated by systematic effects such as the calibration of the jet energy scale. 



3 Point Gib 

In this case the NLSP is neutral and long lived. Almost all of the produced NLSP's exit the 
detector without interacting, so the phenomenology is qualitatively similar to the SUGRA 
models. We first see evidence for new physics via the presence of events with large total 
Et, large IpT-, and isolated leptons. Approximately 0.1% of the NLSP's will decay within 
the detector volume, resulting in a photon that does not point to the interaction vertex. 
The ability to identify such photons and measure the decay point would provide a valuable 
constraint on the lifetime and hence information on C grav . We anticipate that a decay 
probability of 0.1% will be difficult but perhaps possible to detect. Apart from this feature, 
this case is similar to one of the SUGRA cases (Point 4) studied earlier 0,0- In that 
case, there was more structure in the dilepton spectra as significant production of xt an d x\ 
occurs in the decay of a gluino. 



3.1 Effective Mass Analysis at Point Gib 

Events that have at least four jets are selected and the scalar sum of the pt's of the four 
hardest jets and the missing transverse energy fa, formed 

M eS =p T ,l+ PT,2 + PT,3 + PT,4 + $T ■ 

Here the jet p^'s have been ordered such that pt,\ is the transverse momentum of the leading 
jet. Figure |10| shows the distribution in M e ff for events where the following cuts have been 
made. 

• E T > 100 GeV; 

• > 4 jets with p T > 50 GeV and p T ,i > 100 GeV; 

• Transverse sphericity St > 0.2; 

• No n or isolated e with p T > 20 GeV and \rj\ < 2.5; 

• E T > 0.2M eff . 

It can be seen clearly from this figure that there is an excess of events at large M e g over that 
expected in the Standard Model, providing clear evidence for new physics. 
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Figure 10: The effective mass distribution showing LHC Point Gib signal and Standard 
Model backgrounds (hatched histogram) 

3.2 Selection of Dilepton events 

We attempt to extract the decay chain xl ~^ ~^ ^ + ^~Xi as follows. Events are selected 
that have M e g > 1000 GeV and $t > 0.1M e g, and two and only two isolated leptons of 
opposite charge with p? > 20 GeV for electrons, > 5 GeV for muons, and \rj\ < 2.5 for 
both. In order to reduce the combinatorial background we again look at the combination 
e + e~ + /i 4 ~fi~ — e^fi^. Figure [XT] shows the dilepton mass distribution. There is a clear 
endpoint at 



1 - _*i = 105.1 GeV. 



We expect this to have a precision of order 0.1% on the position of this end point. The 
Standard Model background shown on this plot reflects the poor statistics of our sample, 
the actual background is much smaller. The figure also shows a small, but statistically 
significant, peak from Z decays. These are arising from the decay X2 ~* Xi^ which has a 
branching ratio of 9%. The fact that this two body decay is of the same order of magnitude 
as the £ + £~Xi decay is strong evidence that the latter is not arising from the three body decay 
X2 ~^ Xi^ + £~ but from the sequence of two 2-body decays. Additional evidence is provided 
by the shape of the edge. We can be confident therefore that the endpoint is measuring the 
combination of masses in the above equation. 

3.3 Extraction of gluino 

We now select events from the previous sample that have at least two jets with p? > 125 
GeV and require that the dilepton mass is below 105 GeV. The invariant mass of the £ + £~jj 
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Figure 11: The dilepton mass distribution at Point Gib. The Standard Model background 
(hatched histogram) is very small. The fluctuations in our event sample for the combination 
shown i.e. e + e~ + — e ± fi T are therefore larger than one would expect in 10 fb -1 . 



system is shown in Figure [12| where again to reduce combinatorial background we look at 



the combination e + e~ + /x 4 fi~ — e ± /i =F : . All combinations of jets with px > 25 GeV are 
shown. The figure shows a broad peak but no clear structure. The events in this plot 
are dominated by decay g — > qqx® — ► qq£ + £~Xi- However the large combinatorial problem 
prevents a kinematic edge from being visible. In order to estimate the sensitivity to the 
gluino mass, we generate another set of events which differ only in that the gluino mass 
has been increased to 800 GeV. The distribution from this sample is shown as the dashed 



histogram in Figure [12]. The total event rate for this sample has been increased by a factor of 
1.24, to compensate for the smaller production rate and to facilitate a comparison of shapes 
of the distributions. It is clear from this figure that these two curves could be distinguished 
and that a constraint on the gluino mass obtained. The different event rates is not directly 
usable as a mass constraint as the branching ratios are not known. 

The signal can be improved somewhat if we restrict the jets to be those that arise from 



bottom quark jets. The resulting distribution is shown in Figure [13] which is identical to that 
for Figure [T^ except that the jets are required to be tagged as b-jets. The mass distribution 
has a clear structure which reflects the kinematic end point of the dijet-dilepton system at 
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M e is the value of the dilepton edge defined above (the b quark mass is neglected in this 
result). This dilepton-dijet endpoint corresponds to 629 GeV. The dashed line on the figure 
shows the result when the gluino mass is increased to 800 GeV; the kinematic end point is 
now at 673 GeV. The event rates in this plot are rather low and to obtain precision using 
this method 100 fb _1 of integrated luminosity will be needed. 

Following in the spirit of Ref. 0, we attempt to look for kinematic structure in the dijet 
mass distribution from the decay g — > qqx® — > qq£ + £~Xi- We begin by selecting events 
with M cS > 1000 GeV and Ifix > 0.1M eff and four isolated leptons. We retain the event 
if these lepton can be grouped into two opposite-sign, same flavor pairs each of which has 
an invariant mass below 105 GeV. We then select the four jets with the highest transverse 
momenta, and pair them up, selecting the combination that gives the smallest value for the 
sum of the two pair masses. A plot of the mass distribution of a dilepton pair and one of 
the jet pairs is similar, with poorer statistics, to that of Figure there is still no clear 
kinematic feature. Figure [TJ] shows the invariant mass of each of the jet pairs(two entries 
per event). The end-point that one expects to see at rrig — m^o is not visible. However 
the shape of this curve is sensitive to the gluino mass as can be seen by comparing the 
dashed histogram which is the result of the same event selection applied to an event sample 
where the gluino mass is raised to 800 GeV. A detailed study [fj of a case similar to this 
concluded that one could measure the mass difference rrig — rn^o with an uncertainty of order 
3% with an event sample approximately 10 times larger than this one (the cross-section is 
approximately 3 times larger due to the smaller gluino mass of 580 GeV and 30 fb _1 of data 
was assumed). We will conservatively assume that we can determine m g — m^o with an 
uncertainty of 50 GeV, i.e. the difference shown in the figures. 



3.4 Extraction of squarks 

In view of the lack of clear structure in the events of the previous section, it will be very 
difficult to extract the cascade decay q — > gq — > qqqx,2 as was possible in the case of Gla. 
We therefore attempt to extract the decay q — > qW — > qWxi- We first select events with 
M e g > 1000 GeV and Ifi? > 0.1M e g- and only two isolated leptons which can form any of 
the following combinations: e + e + , e~e~, fi~ fi~ or e/i (any charges). In addition we 

require that there be two jets each with px > 450 GeV. Few events pass this selection, but in 
those that do 60% of the directly produced supersymmetric particles are squarks. Figure |15| 
shows the lepton jet invariant mass distribution. In order to assess the sensitivity of this 
distribution to the squark mass, another event sample was produced where the squark masses 
were reduced by 50 GeV. The distribution from this sample is shown as the dashed histogram 
in Figure [1^. The shapes are significantly different; the shifting to the left of the dashed 
histogram is symptomatic of the smaller mass of the squark. A Kolmogorov test applied to 
the shape of the histograms in Figure ^ using finer binning and data samples corresponding 
to 30 fb _1 of integrated luminosity indicates that they are distinct with a probability of 92%. 
Sensitivity to squark masses at the 50 GeV level should therefore be possible. 
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Figure 12: The mass distribution of a dilepton and two jets at Point Gib. The dotted 
histogram corresponds to a different event sample where the gluino mass is changed to 800 
GeV, this distribution has been scaled by a factor of 1.24 to facilitate comparison of the 
shapes. The Standard Model background is shown as the hatched histogram. 
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Figure 13: The mass distribution of a dilepton and two b-jets at Point Gib. The dashed 
histogram corresponds to a different event sample where the gluino mass is changed to 800 
GeV. No standard model background events passed the cuts. 
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Figure 14: The mass distribution of two jets at Point Gib. There are two entries per 
event. The dotted histogram corresponds to a different event sample where the gluino mass 
is changed to 800 GeV, this distribution has been scaled by a factor of 1.18 to facilitate 
comparison of the shapes. No Standard Model events passed the cuts. 




Figure 15: The mass distribution of a lepton+ jet system. The dashed histogram corresponds 
to a different event sample where the squark mass has been reduced by 50 GeV. The Standard 
Model background shown as a hatched histogram is dominated by it production. 
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4 Point G2a 



The total supersymmetry supersymmetry cross-section for both cases G2a and G2b is 23 pb, 
larger than that of cases Gl because the squark masses are considerably lower. Since the 
slepton is the NLSP, almost all supersymmetry events contain a pair of sleptons. Most of 
these sleptons are produced in the decays X2 ~> ^ an d Xi ~> l n the case °f Point G2a, 
the NLSP decays with a very short decay length (52 //) to give an additional lepton. If 
the produced sleptons are selectrons or smuons, the high multiplicity of produced isolated 
leptons will provide both a convenient trigger and the first evidence for new physics. Only 
in the case where both sleptons are staus and decay hadronically, will we have to rely upon 
a jet and missing- E? trigger for the first indication of new physics. 



4.1 Dilepton Distributions 

As the branching ratios for Xi ~* @R ~* G£ + £~ are substantial we can attempt to select this 
decay chain by searching for events with isolated leptons and jets as most of the Xi will De 
produced in the decay of strongly interacting sparticles. Events are selected that have at 
least 4 jets with p? > 25 GeV and \t]\ < 2.5 and at least 4 charged tracks associated with 
each jet (to eliminate jets from hadronic tau decays). An M e g- is formed from the scalar sum 
of the transverse energies of the four jets with the largest px and 

M cS = p T) i + p T)2 + PT,3 + PT,4 + $T ■ 

We then require M eS > 400 GeV and E T > 0.2M eff . 

An additional selection requiring two oppositely charged leptons (either electrons or 
muons) with px > 10 GeV and < 2.5 is made and the dilepton mass distribution formed. 
In order to reduce combinatoric background, we form the combination e + e~ + fi + n~ — e ± yu. =F . 



The mass distribution for this combination is shown in Figure [bj. Two edges are visible at 
52 and 175 GeV. corresponding to 



M 2 - Ml = 52.1 GeV 



and 



M? - Mf = 175.9 GeV . 



The Standard Model background shown on this plot appears not to be negligible. However 
this is somewhat misleading as our background sample does not correspond to the statistical 
fluctuations expected for an integrated luminosity of 10 fb -1 . The number of background 
events can be seen more clearly in Figure [T7| which shows the combination e + e~ + fi + fi~ + 
e ± fi T . It can be seen from this figure that this combination has considerable combinatorial 
background in the signal events; the edge at 175 GeV is less clear. From this plot one can 
estimate that, in the mass range 60 to 170 GeV, the true Standard Model background in 



Figure [lq is ± 60; this fluctuation is insignificant compared with the ~ 4000 signal events 
in the same region. 

It is clear that the statistical error on the precision the measurement of the position of 
the lower edge will be small and the systematic errors will dominate. The higher edge has 
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Figure 16: Mm distribution for the flavor subtracted combination e + e~ + ji + n~ — e ± ji T of 
events having two isolated leptons arising at Point G2a. The Standard Model background 
shown does not represent the statistical fluctuations expected in 10 fb _1 of integrated lumi- 
nosity (see text). 
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Figure 17: M u distribution for all oppositely charged dilepton pairs e + e~ + fi + fi~ + e ± ji T 
in events having two isolated leptons arising at Point G2a. The background, mostly from it 
production, is shown as the hatched histogram. 
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Figure 18: As Figure [16] except that the distribution is shown on a logarithmic scale and 
the fits to the edges described in the text are shown. No background is shown on this figure. 



much poorer statistics, so we have used a fit to estimate how well the endpoint might be 
measured. Before any cuts, the signal distribution in Mu has the form 

dN 



dM, 



a 



A(l + M«)fl(M max - Me, 



where M max is the position of the edge. We add to this a background B taken to be constant 
in the vicinity of the edge and then smear with a Gaussian of width om to be determined 
by the fit. We then use the PAW version of MINUIT [|l6| to fit the data set of Figure |l|, 
whose statistical fluctuations represent approximately 10 fb _1 of integrated luminosity, over 
the ranges 20 GeV < M u < 80 GeV and 150 GeV < M u < 250 GeV. The best fits are shown 



in Figure 18, and the fits give 



52.2661^! GeV 
175.46±°;i GeV 



0"M 

cm 



0.86 GeV 
6.5 GeV 



The errors on M max are determined by MINOS JT6[. We can expect systematic uncertainties 
of order 0.1% on these measurements. We therefore expect that the precision of the upper 
edge will be limited by statistics even for 100 fb _1 of data. For the purposes of parameter 
fitting below, we will assume that the errors are 70 and 270 MeV for 10 fb _1 and 50 and 
180 MeV for 100 fb" 1 respectively. 



4.2 Detection of q R 



Here we use the decay chain — > q\\ — > qi%£ T — > qG£ + £ . The same event selection 
as in Section [O is used with the addition of the requirement that the dilepton pair have 
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Figure 19: M^j mass spectrum at Point G2a for events that contain at least two isolated 
leptons and four jets. The linear fit shown as a dashed line is over the range 390 to 590 GeV. 
The Standard Model background is shown as the hatched histogram. 



mass Mu < 52 GeV and transverse momentum pr(^) > 75 GeV to enhance the probability 
that the leptons come from the same decay chain. The two jets with the largest transverse 
momentum were then selected; each is required to have pt > 50 GeV and and to contain 
least 4 charged tracks. We show the invariant mass of the dilepton and one of these two 
jets in Figure [19] and of one the leptons and a jet in (£ ± j) distributions are show 

in Figure In each case only the smaller of the various mass combinations is plotted 
and again we shown the subtracted distributions e + e~ + n + n~ — e ± / u T which have a cleaner 
structure. 

The £ + £~q distribution has an expected endpoint at 



M max 

1V1 Uq 



Ml -Ml = 640.12 GeV 

qr e R 



while the t^q distribution has an expected endpoint at 



Ml 

Qr 



M% 



M~ 

Ir 

M 2 



289.16 GeV. 



The plots show a linear fit below the end points. These fits extrapolate to end points slightly 
larger than the actual values. The errors on the precision of these values will be dominated by 
the systematic uncertainties in the jet energy scales estimated to be of order 1%. We expect 
that the systematic uncertainty in the measurement of the ratio of these two end-points will 
be less than this. 
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Figure 20: My mass spectrum at Point G2a. The fits shown as dashed lines correspond to 
two separate linear fits over the ranges is over the ranges 150 to 280 GeV and 305 to 400 
GeV. The former extrapolates to an endpoint at 315 GeV. The Standard Model background 
is shown as the hatched histogram. 



We can now solve for the £r, Xi> and masses in terms of the measured endpoints: 

M m 



M; 



Mj lq - M% - Ml 



Mo, 



Mi 



Ml + Mf R 
Mi + Ml 



Note that this method for extracting masses requires only the existence of the decay chain; 
the underlying model is not used in the analysis. Of course the interpretation of the masses 
as those of (Jr, £r, etc., is model dependent. 



4.3 Detection of \i decays 

Approximately 50% of the decays of q~L occur to a Xi ■ There is a decay chain starting from 
Xi proceeding through a slepton and a Xi that gives a final state with three isolated leptons 
viz. 



Xi 



v 



X>1 



u£ + £'£+ 



The combined branching ratio is 29%. We begin with the event selection of the previous 
section, requiring that there be at least three isolated leptons in addition to the jets. Events 
are then required to have at least one opposite sign, same flavor pair with an invariant mass 
in the range 40 GeV < M« < 52 GeV so that they are likely to have come from a Xi 
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Figure 21: mass spectrum at Point G2a. The Standard Model background (hatched 

histogram) is very small. 



decay. Any other pairs of leptons of same flavor and opposite sign with Mu < 175 GeV are 
discarded as they are likely to come from x\ decay. The selected dilepton is then combined 



with any other remaining lepton and the invariant mass of the trio is shown in Figure [21]. If 
all three selected leptons arose from the decay of Xi > this distribution would have a linear 
vanishing at 

M 2 ± - M? + M% - M~ = 85.75 GeV 

There is considerable background in this plot. Nevertheless there is clear evidence of struc- 
ture. While it may not prove possible to extract a precision measurement from this distri- 
bution, it provides evidence for the existence of Xi an d a strong consistency check of the 
model. 



5 Point G2b 

In this case the NLSP is the f\, and it has a long lifetime, cr pa 1 km. The decay £ — > f\rl 
is not kinematically allowed, so the cr and fin are also long-lived. Each event will contain 
two of these quasi-stable particles which will appear in the detector as a pair of slow muons. 
These will provide a trigger as the mean velocity is quite large (see below) and the first 
evidence for new physics in this case. 

A small fraction of the sleptons will decay within the detector; significant energy is 
released in the decay so the result is a track which ends somewhere and a decay product 
that points back to this end point. If the sparticle in question is a selectron, the resulting 
electron could be pointed back using the electromagnetic calorimeter in combination with 
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Figure 22: Generated slepton velocity distribution ft at Point G2b. The dotted curve is for 
sleptons with r] < 1. 



information from the central tracker. If the sparticle is a smuon, the resulting muon would 
have to be pointed using the remainder of the tracking volume and the outer muon system. 
If it is a stau, the resulting hadronic decay of the tau might be pointed using a combination 
of the central tracker and the electromagnetic calorimeter.^ It is worth emphasizing the 
importance of measuring this decay length: it is the only way to obtain information on the 
gravitino couplings and the fundamental scale of supersymmetry breaking. 



5.1 Effective mass analysis at Point G2b 

The events can be triggered using the quasi-stable particles that will appear as muons. The 
velocity distribution of these particles is shown in Figure from which it can be seen that 
the mean velocity is greater than 0.9c. Hence many of these should pass the ATLAS level-1 
muon trigger [ [nfl . 

The events also have a large amount of missing Et as measured by the calorimeter. The 



distribution is shown in Figure 23, which has a mean value of 315 GeV. If the measured 



momenta of the sleptons is included, the missing energy is much smaller as can be seen 



from the dotted curve in Figure ^3[ The true missing Et is larger than standard model 
backgrounds due to the larger number of taus and heavy flavors in the SUSY sample. The 
Standard Model background shown on this plot is controlled by the requirement that it 
contain at least two muons. This calorimetric missing Et could also be used as a trigger. 

We begin the analysis by using the effective mass distribution found by selecting events 
that have at least four jets and forming scalar sum of the p^'s of the four hardest jets and 



§To our knowledge no detailed study of the smuon and stau cases has been done. 
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Figure 23: Calorimetric missing transverse energy fop at Point G2b. The dotted curve shows 
the true Ifa including sleptons. The Standard Model background of dimuon events is shown 
as the hatched histogram. 

the missing transverse energy 

M cS = pr,l + PT,2 + PT,3 + PT,4 + $T ■ 

Here the jet pr's have been ordered such that pr,i is the transverse momentum of the leading 
jet. There is no requirement that the jets or missing Et be large enough to provide a trigger; 
we assume that the events are triggered by the muon system. Note that the Standard Model 
background shown on this plot is required to have two muons in it and is suppressed as a 
result. The NLSP's are ignored when making this effective mass variable. The distribution 
shown in Figure has a mean value of 1004 GeV characteristic of the masses of the strongly 
interacting sparticles which dominate the production. 

There is a peak in the effective mass distribution at zero, corresponding to the production 
of events which have little hadronic activity and Ifix- It is due to the direct production of 
NSLP via such processes as qq — > £ + £~ and the direct production of Xi an d Xi which then 
decay to the NLSP. This is shown in Figure |25] which shows the ID codes for the produced 
primary sparticles and demonstrates that they are mainly sleptons and gauginos. We will 
discuss these events in more detail below. 

5.2 Slepton mass determination 

The sleptons (NLSP) are dominantly produced at the end of decay chains and consequently 
the majority of them are fast. The velocity distribution is shown in Figure ^ from which 
the mean velocity can be determined. Time of flight measurements in the muon detector 
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Figure 24: Effective mass distribution not including sleptons at Point G2b. The Standard 
Model background of dimuon events is shown as the hatched histogram. 
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Figure 25: The particle ID codes for the primary particles of events in Figure ^4| which have 
M eS < 100 GeV: left sleptons (31-36), gauginos (30, 39 and 40) and right sleptons (52, 54, 
and 56). 



system can be used to determine this velocity. When this is combined with a measurement 
of momentum in the same system, the mass can be obtained. 

The muon chambers in the ATLAS detector can provide a time-of-flight resolution of 
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Figure 26: Reconstructed slepton masses for time delays 10 ns < At < 50 ns relative to a 
(3 = 1 particle. 



about Ins. For each slepton with \r)\ < 2.5 the time delay relative to (3 = 1 to the outer 
layer of the muon system, taken to be a cylinder with a radius of 10 m and a half-length 
of 20 m, is calculated using the generated momentum and smeared with a Gaussian 1 ns 
resolution. The smeared time delay At and measured momentum are then use to calculate a 
mass. The resulting mass distribution is shown in Figure 126] for sleptons with 10 ns < At < 
50 ns. Raising the lower limit on At improves the mass resolution but reduces the efficiency 
The resolution is never good enough to resolve the f\ and £r masses of 101.35 GeV and 
102.67 GeV. The upper limit on At is somewhat arbitrary; it reflects practical concerns and 
also eliminates sleptons with very low (3 that lose most of their energy in the calorimeter. 
The average of the generated distribution, 102.2 GeV, agrees well with the fitted mean value 
in Figure [26|. It is important to note that this method will provide an mass measurement of 
an average over the fi, (1r and cr masses as this analysis cannot distinguish slepton flavors. 



5.3 Reconstruction of x\-> X2 ano ^ X4 

Since the £r are quasi-stable, the decays Xi ~~ > ^R^ can be fully reconstructed. Events 
are selected that have at least three isolated electrons, muons or quasi-stable sleptons with 
1 77 1 < 2.5 and pr > 10 GeV. The two highest px particles among the sleptons and muons 
{i.e. those particles that penetrate to the muon system) are assumed to be sleptons and are 
assigned their measured p and the slepton mass; the rest are assumed to be muons. The 
Standard Model background is already negligible, so there is no need to make a time-of-flight 
cut to identify the sleptons. Sleptons are then combined with electrons or muons (£) of the 
opposite charge and the resulting mass for all I^IF combinations is plotted in Figure |2_7|; we 
have no way of determining the flavor of a slepton. There are two narrow peaks at the Xi 
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Figure 27: £ ± £ T mass distribution at Point G2b. 
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Figure 28: Same as Figure ^ for events in a ±5 GeV window around the Xi mass peak of 



Figure 27 and the lepton is identified as an electron. 



and x\ masses. The rather strange shape of the Xi peak is a consequence of the fact that 
the splitting between the Xi an d £r is small, so the mass is dominated by the rest mass of 
the £r. There is a small peak at 348 GeV due to the decay of Xi- 

The mass measurement of the slepton itself can now be refined and the smuon and 
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Figure 29: The mass distribution of a slepton and a pair of leptons at Point G2b. Events 
are selected so that a lepton and a slepton reconstruct to the Xi mass (M^o ± 5 GeV). 

selectron separated by using the events in the in the Xi mass peak. By restricting the event 
samples to the cases where the lepton is either an electron or muon, this method can be 
used to provide separate samples of cr and p,R. The analysis of the previous section is 
now repeated for the events within ±5 GeV of the Xi peak where the lepton is an electron. 
The resulting distribution, shown in Figure ^ has a mean quite close to the correct cr 
mass. The statistical error on the mass in from a data sample of 10 fb _1 is ~ 100 MeV. 
Therefore it should be possible to distinguish the average slepton mass as determined in the 
previous section from the cr and JIr masses. The actual errors are likely to be dominated 
by systematic effects estimate to be 0.1%. This will be sufficient to constrain the mass of 
the stau which is present in the average from the separate selectron and smuon masses. 




5.4 Extraction of £l 

We can begin with the reconstructed Xi an d combine it with another charged lepton in an 
attempt to detect the decay chain £i — > Xi^ ~~ > ^r^^-- We select events that have at least 
two muons or stable sleptons with p T > 10 GeV and \r)\ < 2.5 The two highest pt objects 
are assigned to be sleptons and the rest are called muons. Combinations of a slepton and 
either muon or electron are formed that have no net charge and the system is tagged as Xi 
if Mi±p = Myfi ± 5 GeV. This Xi candidate is then combined with another charged lepton 



and the mass of the three lepton system is shown in Figure There is a clear peak at the 
slepton mass of 203 GeV. 

Another feature is also present in this plot. The decay chain xf ^ u£ ^ x\ v ^ ~ > l%t T vt 
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Figure 30: XrQ mass distribution at Point G2b. Events have a lepton and a slepton with 
invariant mass in a ±5 GeV window around the Xi m ass peak of Figure |27]. The dashed 
distribution corresponds to cases where the jet is from a 6-quark. The fit is a Gaussian over 
the range 610 to 640 GeV. 



has a kinematic upper bound for the mass of the £%£ T £ system of 



M 2 ± — Mj + M 2 



134 GeV 



The structure below this end-point is clearly visible in Figure However there is a large 
background so a very accurate measurement will be difficult. Nevertheless, this structure 
measures a combination of v and \x an d will provide a powerful constraint on the model. 



5.5 Reconstruction of Squarks 

At this point squarks are considerably lighter than gluinos, reflecting the fact that m- g / rriq ~ 
\/Wl and N$ > 1, as is usually the case when the NLSP is a slepton. Direct production of 
squarks is dominant and the decay q^ — > Xil is large. Events are selected that have a £1 
mass within 5 GeV of the Xi peak in Figure 27. This pair is then combined with any of the 
four highest pr jets in the event. The resulting mass distribution of the jet-x? system shown 
in Figure |30] has a relatively narrow peak somewhat below the average q~R mass of 648 GeV. 
Some shift is expected since the jets are defined using a small cone size, R = 0.4. 

The decay q~i — > X2Q is significant but not dominant. It can be reconstructed by selecting 
events that have a ££ mass within 5 GeV of the X2 peak in Figure 27 and combining the 
X2 momentum with the momentum of any of the four highest Pt jets. The resulting mass 



distribution of the jet-X2 system is shown in Figure 31; it has a somewhat wider peak than 
that of Figure B0 a little below the average mass of 674 GeV. The signal to background 
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Figure 31: x^'Q mass distribution at point G2b. Events have a lepton and a slepton with 
invariant mass in a ±5 GeV window around the The dashed distribution corresponds to 
cases where the jet is from a 6-quark.The fit is a Gaussian over the range 625 to 700 GeV 



30] , a reflection of the smaller branching ratio qi 



ratio is poorer than in Figure 
(~ 25%). 

It is also possible to reconstruct the 6i 2 squarks, The subset of events in Figure |3l] for 
which the jet is tagged as a b jet is shown in Figure [31] as the dashed curve. No correction of 
the 6-jet energy is done. The b squark peak is clearly visible and is at somewhat lower masses. 
The resolution is insufficient to separate the peaks from b\ and b 2 whose average mass is 647 



GeV and separation is 9 GeV. A subset of the events of Figure |3(j where the quark jet is tagged 
as a 6-jet is shown as the dashed line in Figure |3(| No structure is visible in this case. This 
difference is due to the branching ratios; BR{pi — > Xify = 5.9%, BR{b 2 — > x\b) = 3.6%, 
BR(q R -> xiq) = 94%, while BR{b\ -> x 2 b) = 15.6%, BR(b 2 -> % 2 &) = 10.8%, and 
BR(q~i — ► x^q) = 25%. The mixing between the b-squarks ensures that both can decay to 

X2- 

The wider peak of Figure [IT] relative to that of Figure |3(] can now be understood. It 
is due to the presence of a significant number of 6-squarks in the former distribution. The 
mass differences between q~L and b are not large enough for the peaks to separate and the 
result is a broad distribution. 



5.6 Reconstruction of fr Decays 



The decay Xi t 1± t =f is more difficult to reconstruct than x° — > but it can provide 

information on the gaugino content of the xl- A technique similar to that discussed in 
Refs. [18] and [19] can be used. Hadronic r's are selected with visible pr > 20 GeV and 



r] < 2.5. These are identified by taking jets where the number of charged tracks is less than 
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M T ~(GeV) 

Figure 32: f-r mass distribution at Point G2b. Solid curve: Using visible r momentum. 
The dotted curve is obtained by selecting events where the missing E T is aligned with the 
tau direction by A0 < O.ln and adding $ T to r momentum. 



or equal to three. The simplest approach is to combine the visible r momentum with the 



slepton momentum. The resulting mass distributions are shown in Figure |32| and on an 
expanded scale for masses near the Xi mass in Figure [330. These curves do not have true 
peaks because of the missing z/ T , but they do have fairly sharp endpoints at the Xi an d x\ 
masses. 

If the slepton momenta are included in the calculation of $t and there are no other 
neutrinos, then fj, can be used to determine the true r momentum. Only the highest pt t is 
used, and the angle between it and the $t direction is required to be A0 < 0.17T. The visible 
t momentum is then scaled by a factor 1 + J^t/Et^, and the f-r mass is recomputed. This 



gives the dashed curves in Figures |32] and As expected, including Ifi? not only reduces the 
statistics but also worsens the resolution for Xi, since the r is very soft in this case. However, 
it produces a peak near the right position for the While this peak probably does not 
improve the mass resolution, it adds confidence that one is seeing a two-body resonance. 

5.7 Direct production of electro-weak sparticles 

We now return to the events present in the peak at low M eS shown in Figure |24|. As pointed 
out above these events are due to the direct production of gauginos and sleptons. We begin 



with the event sample used in Figure 28 and remake that distribution with the requirement 



^This plot requires at one least slepton, assumed to be the "muon candidate" with the largest p t , to be 
present in the event. A second slepton may be required to facilitate triggering. Since all events have two 
sleptons, the efficiency for this is quite large. This additional requirement would reduce the rate shown in 
Figure slightly. 
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Figure 33: Same as Figure [32] on a finer scale with smaller bins. 
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Figure 34: Same as Figure £7] except that the additional requirement M ef j < 100 GeV has 
been made. 



that M c fj < 100 GeV. This is shown in Figure [34[ This plot has a very strong peak at 
the mass of Xi an d a weak, though still clear one, at the mass of \2- This higher peak is 
suppressed as the lepton from its decay is contributing to M eff and the cut throws away some 
signal. 

We can also repeat the analysis of section IO, with the addition of the cut M eS < 100 
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Figure 35: The mass of the slepton-lepton-lepton system formed by selecting events in the 
mass range 110-120 GeV from Figure |33] and combining this slepton-lepton pair with and 
additional lepton. 



GeV. We can take the events in the Xi peak (in the mass range 110-120 GeV) of Figure |4 
and then combine that reconstructed Xi momentum with an additional charged lepton. The 
mass distribution of the resulting system is shown in Figure [35|. This plot shows a peak at 
204 GeV which corresponds to the decay £l — > Xi^ ~~ * ^r^ t ^-- The fit shown on this plot is a 
linear combination of a exponential, a Gaussian and constant over the mass range 135-250 
GeV. The slepton peak is visible, although there are many fewer events than in Figure |2T3| . 
The kinematic feature from the decay chain xt — > vt — > Xi u ^ £%£ T u£ is still visible. 



6 Determining SUSY parameters 

Once a number of quantities have been measured, we can attempt to determine the particular 
SUSY model and the values of the parameters. The strategy will be to attempt to perform 
a global fit to the model parameters using all of the available data, much as the Standard 
Model is tested using the W and Z masses and the many quantities precisely measured by 
LEP/SLC. Such a fit is beyond the scope of our work, and we adopt a simpler procedure. We 
assume that from measurements of global parameters such as those discussed in section we 
know the approximate scale of the superpartner masses and have some idea that we might 
be in a GMSB model. The object is then to determine the parameters of that model and 
check its consistency. We must therefore determine the parameters A, M m , N$, tan/? and 
sgn/i and C grav . If we know the value of one gaugino and one squark or slepton mass of 
the first two generations measured at the mass scale M m then N 5 and A are determined. 
Since these are measured at the a lower energy scale, the physical masses also depend on 
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M m via the renormalization group scaling and we need one more measurement to constrain 
it's value. Two gaugino masses and one squark or slepton mass of the first two generations 
suffice, tan f3 can be constrained either from the Higgs mass or from the masses of the third 
generation squarks and sleptons. In the case of model G2b, the splitting between and fi 
constrains it. Additional constraints on tan/3 and on sgn /i arise from the Higgs, Xi an d xt 
masses. The only constraint upon C grav arises from the lifetime of the NLSP. In cases Gla 
and G2a and G2b we have precise measurements of the slepton and X2 an< ^ Xi masses, so 
the less precise measurements of the squark and gluino masses are not useful in determining 
the fundamental parameters; they only provide powerful consistency checks. Indeed the case 
G2b has so many observables, that it is enormously over constrained. 

In addition to the measurements presented above we assume that the lightest Higgs 
boson has its mass determined precisely from its decay to gamma-gamma. We will assume 
two values for its error; ±3 GeV, which we estimate is the current systematic limit in the 
theoretical calculations needed to relate it to the model parameters; and ±100 MeV which 
corresponds to the expected experimental precision. 

Our strategy for determining the parameters is as follows. We choose a point randomly in 
parameter space and compute the spectrum. We assign a probability to this point determined 
from how well it agrees with our "measured quantities" using our estimates of the errors on 
those quantities. The process is repeated for many points and the probabilities used to 
determine the central values of the parameters, their errors and their correlations. N5 is 
treated continuous variable for these purposes. 

At Point Gla the measurements discussed above for 10 (30) fb^ 1 where statistical errors 
will still be important, viz. 





( Ml R 

M-o 




M\ = 


j M l - 


Ki = 




ML = 




/M-o \ 2 

I^Lj = 105.1 ±0.10(0.10) GeV , 



• m h o = 109.47 ± 3 GeV, 
imply that 

• A = 90000 ± 1800(1200) GeV, 

• M m = 500000 ± 150000(120000) GeV, 

• tan/3 = 5.0 ±1.5(0.9), 

• N 5 = 1 ±0.012(0.10). 

sgn /i is determined unambiguously. Some improvement will be possible with greater in- 
tegrated luminosity until the systematic limit is reached. This will occur for 100 fb^ 1 of 
integrated luminosity. Assuming that the errors are then ±100 MeV, ±200 MeV, ±100 
MeV, ±200 MeV and ±3 GeV respectively, the uncertainties on the parameters reduce: 
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• A = 90000 ± 630 GeV, 



• M m = 500000 ± 80000 GeV, 

• tan/5 = 5.0 ±0.3, 

• V 5 = 1 ± 0.008. 

If the error on the Higgs mass is reduced to ±100 MeV, the uncertainty on tan j3 reduces to 
±0.1; the other uncertainties are unchanged. The poorer precision on M m reflects the fact 
that it enters only via the renormalization group evolution and therefore that the observed 
masses depend only logarithmically upon it. 

At Point Gib, we have for 10 fb _1 of integrated luminosity 



A2 



\ 




1- !t±H i ,1 - (^1 = H).k1 ±0.10 GeV. 



nig — m^o = 523 ± 50 GeV, 



m h o = 109.47 ± 3 GeV 



The precision on the second of these numbers can be expected to increase with more inte- 
grated luminosity; the others are systematics limited. These are sufficient only to constrain 
the following with any degree of precision 

• AN 5 = 90000 ± 860 GeV, 

• tan/3 = 5.01} J 

This is due to an accident in our choice of parameters. The position of the kinematic endpoint 



of Figure [TT| is insensitive to variations of the slepton mass, when vn~ lR = ^om^o. For 
our choice of parameters these quantities differ by 0.5 GeV! As Xi an d Xo are almost purely 
gaugino, these relations provide constraints only on the gaugino masses, i.e. on the product 
AA5. If we assume that the we are able to constrain the average light squark mass within 50 



GeV of its nominal value as appears to be possible from the discussion surrounding Figure [L5|, 
for 30 fb _1 , we then obtain 

• AN 5 = 90000 ± 620 GeV, 

• A = 90000 ± 8100 GeV, 

• M m < 7 x 10 8 GeV (95% confidence) 

• tan/3 = 5.0l}j 

Here we have restricted M m > A; the model is not sensible if this is not the case. If the error 
on the Higgs mass is reduced to ±100 MeV, the uncertainty on tan/3 reduces to ±0.1; the 
other uncertainties do not reduce. It is not possible to determine sgn fi using the signals that 
we have shown. For example, the set of parameters A = 104500 GeV, M m = 0.239 x 10 9 
GeV, N 5 = 0.872 and sgn/z = —1 is acceptable. The mass of %\ (x?) is increased by 19 (6) 
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GeV and the mass of 6r by 40 GeV, but this case has the same end point in Figure |TT[ An 
independent constraint on the slepton mass or a measurement of the squark mass with a 
precision of order 10 GeV is needed to eliminate this case. 

At Point G2a the measurements discussed above for 10 fb _1 viz. 



m\ - m \ = 52.21 ± 0.07 GeV, 



m 2 _ m 2 = 175M ± Q GeV, 



m\ R - m\ R = 640 ± 7 GeV, 



m~ — m~ 

1R e R 



\ 



M? 



X 



1 - ^/./M'i -M-L = 0.45 ± 0.004 



■ M in. Xl 



• m h o = 106.61 ± 3 GeV, 
imply that 

• A = 30000 ±410 GeV, 

• M m = 250000 ± 44000 GeV, 

• tan/5 = 5.0 ±0.7 

• iV 5 = 3 ± 0.036. 

sgn fi is determined unambiguously. Only small improvements can be expected as the inte- 
grated luminosity is increased above 10 fb _1 . If we reduce the errors on the first two quantities 
to 50 and 180 MeV respectively, their likely systematic limits, as might be achieved with 30 
fb _1 of data we obtain 

• A = 30000 ± 380 GeV, 

• M m = 250000 ± 42000 GeV, 

• tan/3 = 5.0 ±0.7 

• N 5 = 3 ± 0.035. 

If the error on the Higgs mass were reduced to ±100 MeV, the error on tan/3 reduces to 
±0.04. 

At Point G2b the measurements discussed above viz. 

• m~ eR = 102.67 ±0.1 GeV, 

• m~ m = 102.67 ±0.1 GeV, 

• m fl = 101.35 ±0.1 GeV, 

• m,o = 115.18 ±0.1 GeV, 
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• m^o = 203.71 ± 0.2 GeV, and 

• m h o = 106.61 ± 3 GeV, 
imply that 

• A = 30000 ± 175 GeV 

• M m = 250000 ± 22500 

• tan/? = 5.0 ±0.21 

• N 5 = 3 ± 0.014. 

sgn /i is determined unambiguously. These measurements are likely to be at their systematic 
limits with 10 fb _1 of integrated luminosity and further improvement will be difficult. If the 
error on the Higgs mass is reduced to ±100 MeV, the error on tan/3 reduces to ±0.02. 

In case Gib whose global signatures are similar to SUGRA, we must address the issue 
of whether it could be confused with a SUGRA model. We search SUGRA parameter space 
for test of parameters that could be consistent with the following: 



• rrig — m^o = 523 ± 50 GeV, 

• m h0 = 109.47 ± 3 GeV, 

• nig -m~ g = 200 ± 500 GeV. 

We obtain the following solution for the SUGRA parameters 

• m = 100 ± 20 GeV 

• mi/a = 295 ± 6 GeV 

• tan/3 = 4.5 ± 1.1 

• sgn fj, = +1 

• A = 250 ± 200 GeV 

This solution has only a 15% probability. The central value of the light squark masses for 
this solution is 760 GeV. Most of the other masses are similar to those of case Gib with the 
exception of cl which has mass 100 GeV larger. This result illustrates the general difference 
between SUGRA and GMSB models. The mass splitting between the squarks and e# is larger 
in the GMSB case. If we assume that rrig — m g = 200 ± 75 GeV as we might expect from 
the methods of Section || with at least 30 fb -1 of integrated luminosity, then the SUGRA 
solution is eliminated (it has 10 -5 probability) and the ambiguity resolved. 
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7 Discussion and Conclusions 



In this paper we have given examples of how LHC experiments might analyze supersymmetry 
events if SUSY exists and if the pattern of superpartner masses is given by gauge mediated 
models of supersymmetry breaking. We have illustrated the four classes of phenomenology 
to be expected in such models: events with missing energy that are similar to those expected 
in SUGRA models (Point Gib); events with a pair of isolated photons from xi decay (Point 
Gla); events with long lived sleptons (Point G2b); and events with leptons from prompt 
slepton decay (Point G2a). In the first case, we have discussed how measurements can be 
made which enable one to prove that the gauge mediated model and not SUGRA is respon- 
sible for the pattern of masses. In the other cases, detection and measurement is easier than 
in SUGRA. Characteristic features are present such as photons, leptons or stable charged 
particles, that enable the supersymmetry signals to be extracted trivially from Standard 
Model backgrounds. In addition, these features make it possible to identify and use longer 
decay chains. 

We have illustrated a technique whereby the supersymmetry events can be fully recon- 
structed despite the presence of two undetected particles. The technique relies upon there 
being a decay chain of sufficient length that occurs twice in the same event. Each step 
in the chain provides a constraint and sufficient constraints can occur that together with 
the measurement of ft? enables the event to be reconstructed. In such cases the masses of 
the superparticles in the decay chain can be measured directly rather than inferred from 
kinematic distributions. 

In all of the cases discussed in this paper, as in the SUGRA cases discussed previously 
[H, 0, HI, the LHC will be capable of making many precise measurements that will enable 
the underlying model of supersymmetry to be severely constrained should supersymmetric 
particles be observed. The key to all of these analyses is the ability to identify a characteristic 
final state arising from the decay of a sparticle that is copiously produced. The main such 
decays that have been used in the SUGRA and GMSB analyses done so far are: 

• Dileptons from \2 — > Xi^ + ^~ or X2 — ► £r^rXi^ + ^i 

• taus from the decays \2 — > Xi'h +T i or Xi T\T\ —* Xi T+T ~ > which can dominate 
when tan/3 is large; 

• X2 -> Xih -> X°ibb. 

In the case of GMSB or i?-parity breaking models the subsequent decay of Xi can provide 
additional information and constraints. Of course the information extracted in this way 
is only a small fraction of the total available. A complete analysis will involve generating 
large samples of events for many SUSY models and comparing all possible distributions with 
experiment. 
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Appendix 



The details of the full reconstruction used in section ^72] are given here. 

Events are selected to have four leptons and two photons with a unique combination of 
two leptons and one photon coming from each x\ decay. The xHi £r, and Xi masses are 



assumed to be known precisely from the distributions discussed in section |27i], and the G 
mass was assumed to be (essentially) zero. This leads to the following set of equations for 
the 4-momentum p of the gravitino: 

p 2 = 

(p + p y ) 2 = Mfo, 

(p + P-t+Pi,) 2 = Mf, 
(p + P 7 +Pe 2 +Pe 1 ) 2 = M|o. 

This implies that the gravitino momentum p is given in terms of the photon momentum k 
and the lepton momenta I and q by 



2p k 


-2p-k = 


M h 






2p l 


-2p-l = 


M l 




-2k-l = Q, 


2p <?o 


— 2p ■ q = 


M h 




-2(k + l)-q-q 2 = C q , 



where 

Po = \JP- 

It is now clear that they give two linear and one quadratic constraint and hence an additional 
2x2 ambiguity. The solution is straightforward. From the above equations one finds 

p-Dt = E lt 
P-D 2 = E 2 , 



where 



These can be solved to give 
where 



F. 



Di = 2l k - 2k l , 

Ei = —loCk + k Ci , 

D 2 = 2q k - 2k q, 

E 1 = -q C k + k C q . 

Pi = Fi + Gip z , £ = 1,2, 

E\D 2 y — E 2 D ly 
D\ x D 2v — D 2x D\ 



'2y — U2x±J\y 

n _ Di z D 2y — D 2z Di y 

{j-r — ^ ^ — — 



D\ X D 2 y — D 2x D\y 
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E\D2x — E 2 Di x 

D\yD2x — D2yD\ x ' 
DlzD2x — D2 z Dlx 

DiyD 2x — D 2y D lx 

This is then substituted into the first of the original equations to yield a quadratic equation 
for p z , 

H + H lPz + H2p 2 z = 0, 

where 

H = 4k 2 (F 2 + F y 2 + M 2 o)-(C k + 2F x k x + 2F y k y ) 2 , 

H x = 4k 2 (2F x G x + 2F y G y )-2(C k + 2F x k x + 2F y k y )(2G x k x + 2G y k y + 2k z ), 
H 2 = 4k 2 (G 2 x + G 2 y + l)-(2G x k x + 2G y k y + 2k z ) 2 . 

This gives two solutions for the G momentum p for a given assignment of the other momenta 
provided Hf — AH H 2 > and no solution otherwise. 
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